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Abstract 

Recently it has been proposed that a strange logarithmic expres- 
sion for the so-called Barbero-Immirzi parameter, which is one of the 
ingredients that are necessary for loop quantum gravity (LQG) to pre- 
dict the correct black hole entropy, is not a sign of an inconsistency of 
this approach to quantization of general relativity, but is a meaningful 
number that can be independently justified in classical GR. The alter- 
native justification involves the knowledge of the real part of the fre- 
quencies of black hole quasinormal modes whose imaginary part blows 
up. In this paper we present an analytical derivation of the states with 
frequencies approaching a large imaginary number plus ln3/87rGArM; 
this constant has been only known numerically so far. We discuss the 
structure of the quasinormal modes for perturbations of various spin. 
Possible implications of these states for thermal physics of black holes 
and quantum gravity are mentioned and interpreted in a new way. A 
general conjecture about the asymptotic states is stated. 



e-print archive: http://lanl.arXiv.org/abs/gr-qc/0212096 
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1 Introduction 

While there are several open issues about the way quantum dynamics is 
currently handled in loop quantum gravity, the pretty kinematics — or per- 
haps the mathematically elegant attempts to understand de Sitter space in 
this framework (see [1]) — can be separated and could well serve as a useful 
piece of our future understanding of quantum gravity without recourse to a 
background. 

Under some assumptions, loop quantum gravity implies that the black 
hole entropy is proportional to the area of the horizon ^ [2, 3] as well as a 
numerical constant called the Barbero-Immirzi parameter [4, 5] whose value 
must be chosen appropriately so that the result agrees with the Bekenstein- 
Hawking entropy. The frustrating role of the Barbero-Immirzi parameter 
recently changed a bit because of an interesting paper by Dreyer [6]. Build- 
ing on a proposal by Bekenstein and Mukhanov [7] and extending a numer- 
ical observation by Hod [8], Dreyer proposed that the gauge group of LQG 
should be 5*0(3) rather than SU{2). Such a minor modification changes the 
Barbero-Immirzi parameter required to reproduce the Bekcnstein-Hawking 
entropy from 75^/(2) = lii(2)/(7rA/3) to 7^0(3) = ln(3)/(27r-\/2). More impor- 
tantly, this new constant leads to a new value of the minimal positive area 
Aq = 4:Gn In 3. The frequency emitted by a black hole which is changing its 
area by this amount turns out to coincide with another frequency that has 
been only known numerically: the asymptotic real part of the frequencies of 
black hole's quasinormal modes. ^ They will be the focus of our paper. 

Although an immediate reaction is that this agreement is a coincidence, 
we will present some evidence that it does not have to be an accident. For 
example, a concern is that the real part agrees with the conjectured value 
ln(3)THawking Only in the first few decimal places and the following ones 
will show a discrepancy. We will show that this is not the case. The most 
important result of this paper is an analytical proof that the asymptotic 
frequency is indeed equal to ln(3)rHawking for perturbations of even spin. We 
will also present several other arguments that the agreement might indicate 
something important. We will use methods of classical general relativity. 
Every good theory of gravity should probably reproduce these results (for 

^We were informed that this idea of Dreyer was actually first proposed by Kirill Kras- 
nov in an unpublished email correspondence between K. Krasnov, J. Baez, and A. Ashtekar 
in 1999, soon after Hod's paper [8]. The physicists stopped thinking about this proposal 
because of two reasons. The continuous character of the Hawking radiation did not agree 
with their attempts to interpret the frequency of the quasinormal modes. We will men- 
tion this problem in subsection 2.2. The main reason was puzzling numerical results by 
Kokkotas et al. [9] that we will discuss in subsection 3.3. 
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large black holes, the frequencies can be still taken to be much smaller than 
JTipianck, for example). The thermodynamic nature of our answer suggests 
that the result might tell us something nontrivial about quantum gravity. 



2 The quantum mechanical problem 



The possible radial spin-j perturbations of the four-dimensional Schwarz- 
schild background [10] are governed by the following master differential equa- 
tion [11]: 

■^,+v{x)-u.^y=o. (1) 

We will treat this equation as a Schrodinger equation with the (Regge- 
Wheeler) potential V{x): 

Morally speaking, the perturbation contains a factor of the spherical har- 
monic Yifn with the orbital angular momentum equal to / that appears in 
the "centrifugal" term of (2). The numerator (1 — j"^) of the other term is 
usually called a and is characteristic for perturbations of various spins. It 
might be useful to summarize its values for the most important examples of 
j- 



1 

-3 



scalar perturbation J = 

electromagnetic perturbation j = 1 (3) 
gravitational perturbation j = 2 



The "tortoise" coordinate x is related to r by 



d ( \\ d 
x = ln(r-l) + r, ^ = 1 - - W". (4) 



dx \ r J dr 

We chose units with unit radius of the black hole's horizon: 2GnM = 1. 
Therefore the difference (1 — 1/r) is simply the usual Schwarzschild warp 
factor. Also, in the equivalent quantum mechanical problem (1) the doubled 
mass equals 2m = 1. 



We want to study the equation (1) on the interval 
X G (— oo, -|-oo) i.e. rG(l,oo). 



(5) 
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Because the potential is mostly positive (for j < 1 it is strictly positive), there 
are no discrete normalizable bound states. Nevertheless it makes mathemat- 
ical sense to look for discrete quasinormal states, analogous to quasistation- 
ary states in quantum mechanics whose frequency is allowed to be complex. 
These states are required to have purely outgoing boundary conditions both 
at the horizon (r = 1) and in the asymptotic region (r = oo): 

(^(ar) ~ c±e^^'^^ for x ^ ±oo. (6) 

Such a constraint is as strong as the requirement of normalizability and 
singles out discrete solutions u. It is however very hard to look for the right 
wave functions numerically because (6) requires the exponentially decreasing 
component of the wave function to be absent asymptotically. Separating the 
exponentially small contribution from the rest is an extremely difficult task, 
especially for large decay rates. 

The time dependence of the wave function in our quantum mechanical 
model is then e*'^*, and because we want to study exponentially decaying^ 
modes, the imaginary part of uj will always be taken positive. This also 
implies that while the sohitions can oscillate, they must also exponentially 
increase with There is a reflection symmetry 

LO ^ —LO* (7) 

which changes the sign of Re(a;). Although the usual conventions describe 
the wave e"*"^^ as outgoing at a; ~ +oo for Re(aj) < 0, we will use the 
reflection symmetry and assume, without loss of generality, Re(a;) > 0. On 
the other hand, the assumption Im(ti;) > is physical and important because 
our final result will originate, in a sense, from positive integers n such that 
(n + 2iuj) 1. 

While the quasinormal modes with small enough lm(ui) represent the 
actual damped vibrations of the black hole — those that bring it quickly to 
the perfect spherical shape, the physical interpretation of the very unstable, 
high-overtone states with |Im(aj)| > Re(aj) is more problematic. Despite 
these conceptual problems, the existence of such formal solutions implies 
the presence of poles in the transmission amplitude (or equivalently, the 
one-particle Green's function). Such poles can be important just like in 
many other situations and we will investigate some possible consequences at 
the end of this paper. 

•^The relation between the energy in the quantum mechanical model and the spacetime 
energy is indirect: note that the former equals while the latter is equal to w. The 
spacetime equations of motion are second order equations. Therefore one should not 
worry about our nonstandard sign in e"^* that was chosen to agree with literature. 
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The following two subsections review the existing literature on the quasi- 
normal modes and basics of LQG. The reader interested in our analytical 
computation is advised to continue with section 3. 



2.1 Existing calculations 

Numerically it has been found by Leaver [12], NoUert [13], Andersson [14] 
and others that the n**^ quasinormal mode of a gravitational j = 2 (or scalar 
j = 0) perturbation has frequency 

i(n-l/2) ln(3) ^, _i/2, 
0Jn= ^ ^ ' ' +^ + 0{n ^/^), n^oo. (8) 

This asymptotic behavior is independent of the orbital angular momentum I 
although the detailed values are irregular and ^-dependent for finite Im(a;). 
The asymptotic value (8) is valid for j = 2 as well as j = 0, for example; we 
will argue that Re(ti;) must asymptotically vanish for odd and half-integer 
values of j. The number of quasinormal modes is infinite. This fact was 
proved ten years ago [15] by methods that did not allow the authors to say 
anything about the frequencies. The real part in (8) equals 

Re{ojn) «i 0.087424. (9) 

47r 

This constant has been determined numerically with the indicated accuracy. 
Such a frequency is related to a natural value of the minimal area quantum 
41n(3)Gjv, one of the values 41n(A;)GAr favored by Bekenstein and Mukhanov 
[7] because of their ability to describe the black hole entropy from area- 
degenerate states.^ Four years ago, Hod [8] became the first one to notice 
that the constant (9) can be written in terms of ln(3) and proposed a heuristic 
picture trying to explain this fact. As far as we know, our paper contains 
the first analytical proof that this number is indeed ln(3)/(47r). 

Guinn et al. [18] presented two similar and very accurate higher-order 

WKB calculations. However, an undetermined subtlety makes their results 
invalid because they obtain instead of ln(3)/(47r) for the real part. The 
very recent paper [19] can be viewed as a refined version of [18]; it leads to 
results that are compatible with the present paper. 

^The idea of an equally spaced black hole spectrum goes back to 1974 when it was 
first proposed by Bekenstein. While Bekenstein and Mukhanov prefer the value k = 2 
[16] which also agrees with Wheeler's idea It from Bit, Hod found the value fc = 3 more 
reasonable. Also the authors of [17] propose some arguments in favor of fc = 3. Although 
we disagree with the proposal of equally spaced black hole spectrum, we are grateful to 
J. Bekenstein for correspondence about the early history of black hole thermodynamics. 
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2.2 Loop quantum gravity and Dreyer's observation 

Ashtekar proposed new variables for canonical quantization of Einstein's 
theory. See for example [20] for an efficient review of Loop Quantum Gravity 
or [21] for a more extensive one. The spatial three-dimensional components 
g"-^ of the metric at t = are expressed in terms of the inverse densitized 
dreibeins 

3 

(det53x3)5"' = E^"^' (10) 

i=l 

The index i = 1, 2, 3 is an adjoint index of some 50(3) (or formerly SU{2)) 
gauge theory and the components i?? themselves are interpreted as the dual 
variables to the gauge field: 

Ef (X) = -8"7G»;^ (11) 

Note that the case of three spatial dimensions is special because the adjoint 
and vector representations of 50(3) have the same dimension. The Hilbert 
space is a priori the space of functional of Al^{x). However it is possible 
to consider the subspace of functions of S'0(3)-valued open Wilson lines 
defined on links of a "spin network" ; this is known as loop transform. The 
union of such subspaces for all possible spin networks is dense in the original 
Hilbert space. The functions on the Cartesian product of many copies of the 
50(3) group manifold can be expanded into spherical harmonics. The open 
Wilson lines, i.e. the links, then carry a spin Jj associated with the spherical 
harmonic, "coloring" the spin network. Colored spin networks then form an 
orthogonal basis of the Hilbert space. 



The most interesting result of this approach to quantization of gravity is 
the following area quantization law [22, 23] that results directly from (10) 
and (11): 

A = SttGnI J2 VMJ^+^ (12) 

i 

Here Jj = 0, 1, 2, . . . are spins defined on the links i of a spin network (not to 
be confused with j introduced in (2)) that intersect a given two-dimensional 
sheet in the coordinate space. The physical area A of the sheet depends 
on the dynamical metric variables. According to LQG, the physical area 
is "concentrated" in the intersections with the spin network. Recall that 
colored spin networks form a basis of the Hilbert space of an 50(3) gauge 
theory.^ The constant 7 denotes the so-called Barbero-Immirzi parameter 

''Originally, Ashtekar and Baez worked with an SU{2) theory instead of SO{3). They 
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[4, 5], a pure number that can be interpreted as the finite renormahzation 
of Gn between the Planck scale and low energies. A specific value of 7 is 
required for LQG to reproduce the Bekenstein-Hawking entropy as we will 
recall soon. 

The area (12) becomes a quasicontinuous variable if it is large because of 
the contributions with J > 1: the ratios between numbers Ji{Ji + 1) for 
different Ji are irrational. Nevertheless the typical microstates of a horizon 
are dominated by links with the minimal possible value of J. If the links 
with higher values of J were absent altogether, the spectrum of the area 
operator (12) would become equally spaced. The step would be equal to 
4Gn ln(2J^™°) + 1) for the entropy to come out correctly [7]. If the relation 
between the area and the mass is preserved (which is however hardly the 
case microscopically), the links with J > 1 are essential to keep the energy 
spectrum as well as the spectrum of Hawking radiation continuous. Further- 
more the emission of a typical Hawking quantum should rearrange a large 
number of links — including links with J > — if LQG is to be compat- 

ible with Hawking's semiclassical results, otherwise the spectrum would be 
peaked around discrete frequencies. This is one of the problems that was 
mentioned in footnote in the introduction. We believe that the assumption 
that the energy spectrum is discrete with spacing ln(3)/(87rGArM) was in- 
compatible with the spectrum of Hawking radiation whose energies are of 
the same order but continuous. The correct interpretation is probably more 
subtle. If we reject Hawking's prediction of highly thermal radiation, then 
we also have no more reasons to believe the calculation of black hole entropy 
that is an implication of the temperature calculation; at least the numerical 
coefficient would be unjustified. We find it hardly acceptable to reject Hawk- 
ing's semiclassical calculations, in part because they have been confirmed by 
many developments in String Theory. 

The need for a macroscopic rearrangement of many links (including those 
with J > j(™™)) in LQG to explain a single Hawking particle should not be 
too surprising because LQG needs a similar miraculous interplay between 
a large number of links to generate the flat space itself. One should also 
note that the problem is equally serious regardless of the value of and 
therefore we should not discard Dreyer's proposal to switch from SU (2) to 
SO{3) on these grounds. 

were motivated by the possibility of incorporating fermions into a future version of loop 

quantum gravity. Wc believe that there is no reason why a theory that is meant to rewrite 
gravitational physics should "know" about our desire to incorporate fermions. Therefore, 
SO(3) is more natural than SU{2). The SU{2) gauge theory could lead to similar physics 
but Ji could also be half-integers in this case. The resulting SU{2) Barbero-Immirzi 
parameter ln(2) / (ttVs) has not been supported by any numerical coincidences. 
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While loop quantum gravity does not yet offer a well-defined framework 
to explain the key question why black hole entropy has no volume extensive 
contributions coming from the interior, there exist arguments that LQG im- 
plies the correct entropy of the horizons of all possible black holes i.e. the 
entropy proportional to their area (up to the universal constant mentioned 
a moment ago). See [2] and [3] as well as references therein for more infor- 
mation. Instead of the sophisticated Chern-Simons computations initiated 
by Krasnov, let us review a simpler version of the argument. The observ- 
ables J-^^ living on the links intersecting a given area become undetermined 
and physical if the area happens to be a horizon. The mixed state of the 
black hole is dominated by microstates with the minimal allowed value of 
j(min) _ occupying most links (for the SU (2) gauge group, we would have 
j(mm) _ 1/2). Because of (12), the number of such links is 

A A 

iVlinks = = . (13) 

Each link carries (2j(™'")-|-l) possible z-components of the spin and therefore 
the number of microstates equals the corresponding power of (2j(™"^) + l); 
namely 



-^microstates — GXp 



ln(2j(°^^'^) + 1)- ^ 



87rG7v7VJ(™'i)(j(>^i°) + 1) 
This is equal to the exponential of the black hole entropy A/AGn for 



(14) 



ln(2 J(°^^°) + 1) 



27ri/j(™''i)(j(™°) -I- 1) 



j(min)^;^ 



Such a value of 7 may seem very unnatural. However Dreyer [6] has recently 
suggested an independent argument for this particular value of 7. Because 
we disagree with several details of his argument, we will present a modified 
version of it. 

Assume that a new link with = 1 is absorbed by a black hole 

horizon (or it is created there). Its area therefore increases by Aq, which is 
according to (12) and (15) equal to 

A^ = = 41n(3)Gjv. (16) 

Because A = AttR'^ = 16tt{G]\[M)^, this variation can be written as 
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It follows there should be a preferred frequency (energy) 

unn^ = ln(3)/(87rG^^M) (18) 

associated with the emission of the most generic link. While black holes emit 
Hawking radiation with continuous frequencies, the special frequency tJunk 
turns out to be exactly equal to the asymptotic real part of the quasinormal 
frequencies in (8) (divided by 2GnM = 1). We believe that the special 
frequency should not be interpreted as the energy of a physical Hawking 
quantum which is continuous. 

Because Dreyer was forced to switch from SU{2) to 50(3) and because 
the justification of ln(3) on one side was only numerical, it is clear that his ar- 
gument is nontrivial. Although we do not think that the imaginary part can 
be ignored in Dreyer's mechanistic arguments involving the Bohr's frequency 
and one can criticize other issues, such a coincidence nevertheless deserves 
an explanation. The first step is to extend the numerical evidence for the 
ln(3) form of this asymptotic frequency into an exact analytical derivation. 
We will realize this step in the next section. 

3 The asymptotic expansion for large frequencies 

The solution of (1) can be expanded in the following way (see for example 
[13]): 

*M=(^)'"e-'-*'-'f;<.„(^)". (19) 
The prefactor is meant to satisfy the boundary conditions (6) : 

• e~*'^(''^^) is necessary for the correct leading behavior at x — oo i.e. 

r ^ oo; 

• (r — 1)"^ fixes the behavior at x — — oo i.e. r — 1+; 

• r""^^ repairs the subleading behavior at r — ^ oo coming from the 
logarithmic term in (4). 

The power series (19) converges for 1/2 < r < oo and also for r = oo if the 
boundary conditions at r = oo are preserved. The equation (1) is equivalent 
to the recursion relation 



CQ{n,(jj)an + ci{n,uj)an-i + C2(n,u;)a„_2 = 



(20) 
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where the coefficients Ck{n, uj) can be extracted from the equation and rewrit- 
ten in the following convenient way: 

co{n,uj) = +n{n + 2iLo) (21) 
ci(n,a;) = -2{n + 2iuj - l/2f - l{l + 1) + f - 1/2 (22) 
C2{n,uj) = +{n + 2iuj-lf-f (23) 

Note that except for the factor n in Co(ra), the coefficients Cfc(n, cu) depend on 
n,uj only through their combination (n + 2iu;). The initial conditions for the 
recursion relation guarantee that the boundary conditions at r = 1 will be 
preserved. They are Oq = 1 (or any nonzero constant) and a_i = (which 
also implies a_„ = for all positive integers n as a result of (20)). We also 
define 

Rn = -^. (24) 

«n-l 

The minus sign was chosen to agree with [13]. The relation (20) then implies 
ci(n, UJ) - co(n, ujRn = — (25) 

or, equivalently, 

= f"--^ (26) 

Rn can be consequently written as an infinite continued fraction. The bound- 
ary condition at r = 1, i.e. a_i = 0, take the form 

Ro = oo i.e. ci(l,a;) - co(l,a;)i?i = 0. (27) 

The convergence of (19) at r = oo implies a particular asymptotic form of 
Rn for large n (and < 1 for n oo) and the boundary conditions 

require (27). The equation (27), written in terms of continued fractions, is 
then a condition that oj must satisfy for a quasinormal mode to exist: 

= ci{l,iu)-coil,Lu) ^^^^ J—. (28) 

Ci (2, U) - Co(2, ^) ,^(3,,)_,;(3,J)^miin... 



3.1 Going to the limit 

Let us now consider uj with a huge positive imaginary part. Because Rn — ^ 
— 1 as n ^ oo [13] (in other words, an/an-i — 1) and because for large 

the changes of i?„ slow down, it is an extremely good approximation to 
assume that Rn changes adiabatically as long as n, and \n + 2iu\ 3> 1: 

-^ = l + 0(|a;|-n (29) 
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Such an approximation is good for Re(n + 2iuj) > (where we calculate i?„ 
recursively from R^o = —1) as well as Re(ra + 2iLo) < (where we start the 
computation with Rq = oo). Note that a similar fact holds for the coefficients 
of the Taylor expansion of the exponential. By inserting Rn-i = Rn into 
(25) we obtain a quadratic equation^ with the following solutions (we keep 
only the leading terms for n ^ \2u\ — > oo): 

Rt = + + + 0(|a;|-V2). (30) 

This approximation for the continued fractions defined in (26) has been 

successfully checked with Mathematica. Irregularities in i?„ would imply 
violations of the boundary conditions although we do not claim to have a 
totally rigorous proof. Rn should satisfy (30) for one of the signs; it is a 
necessary condition and a more detailed discussion whether the condition is 
also sufficient is needed. For Re(n + 2iuj) < the sign is determined by the 
condition Ri is small; two terms in (30) must approximately cancel. The 
sign for Re(n + 2iu;) > follows from the requirement \Rn\ < 1 for n — > oo 
which is necessary for the convergence of (19) for r = oo. 

For very large |a;|, we can choose integer M such that 

l<M<|a;|. (31) 

For such intermediate values of M the equation (30) can still be trusted 
when we want^ to determine R[-2iu;]±M'-i a-iid only the second term of (30) is 
relevant in the limit (31). This term makes the numbers R[-2ioj]-\-x scale like 
:^i\fx I \/ —2iuj for moderately small x such that 1 ^ x <C |w| while the first 
term in (30) scales like xjuj and is sublcading. Note that neglecting this first 
term is equivalent to neglecting ci(n,u;) in the original quadratic equation 
(25) ; this term will turn out to be irrelevant for all values of n in the large 
|a;| limit, with the possible exception of some purely imaginary frequencies 
where co(n) or c^iji) vanish. The ratio of -R[-2iw]±M can be calculated from 
(30) as the ratio of for x = ztM: 

§=^^^^ = ±z + 0(M-V2) (32) 

^Matthew Schwartz independently realized that this quadratic equation gives a good 
estimate of the remainders. In fact, he even showed that the solution (30) of this quadratic 
equation coincides with the fully resummed expression by NoUert [13] in the scaling limit 

n ~ \2iuj\ oo. 

''The symbol for the integer part [— 2itt;] should be actually understood as an arbitrary 
integer that differs from —2iu! by a number much smaller than M; we will assume that it 
is even. 
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Because the assumptions that led to (30) break down for \n + 2iuj\ ~ 1 
where the coefficients co(ra), ci(n), C2(n) strongly depend on n and further- 
more contain terms of order 1 (or j^) that now cannot be neglected/ the 
quantities -R[_2ia)]+M ^^'^ R[-2iu!]-M must be related through the original 
continued fraction (instead of the adiabatic approximation which breaks 
down in this region). In the next paragraph we will be able to calculate 
the continued fraction exactly in the |ti;| oo limit. The continued fraction 
must predict the same ratio (32) as the adiabatic argument because these 
two solutions must be eventually "connected." This agreement will impose 
a non-trivial constraint on lo. 



The scaling of R[-2iui]+x as l/-y|tJ| (see (30)) also implies that ci(n) in 
the denominator of (26) (which is of order one) is negligible if compared to 
the other term (which scales like l/R[-2iu)]+x or \/M")- Keeping M fixed and 
scaling — > oo, it is clear that the effect of ci(n) in (26) goes to zero. The 
only possible exception is the case when co(n) and/or C2(n) become almost 
zero for some n; this exception cannot occur for nonzero Re(a;) but we will 
study it in the next subsection. 



Because ci(n) does not affect the asymptotic frequencies, the orbital 
angular momentum I is irrelevant, too. Therefore the asymptotic frequencies 
will depend on j only (which is contained in the definition (23) of C2(n)). The 
factor n in (21) can be replaced by the constant [— 2iu;] as the frequency goes 
to infinity if we want to study a relatively small neighborhood of n ~ [— 2ia;] 
only. Moreover, the continued fraction simplifies into an ordinary fraction. 
Inserting (26) recursively into itself 2M times tells us that 

_ ^ f C2{[-2iuj] - M + 2k - l)co{[-2iuj] - M + 2k) \ ^ 
ii[-2M-M - n (^co([-2M -M + 2k- 1)C2([-2H -M + 2k) ) ^[-^-l+M 

(33) 

The dependence of all coefficients Cfc(n) on uj was suppressed. Now we can 
combine (32) and (33) to eliminate the remainders we will formulate 
the requirement that the generic adiabatic solution (30) that is valid almost 
everywhere can be "patched" with the exact solution for n ~ [— 2za;] (33) 
that has also been simplified by the large |a;| limit. 

Fortunately, the products of co(n) and C2(n) can be expressed in the 
language of the V function (partially because — p = {x + j){x — j); this is 
a deeper reason why we write o" as 1 — j"^). Since C2(n) is bilinear, the four 



'^Except for the discrete character of n in the present case, the breakdown of the generic 
result (30) is analogous to the way how the simple Ansatz for the wave function in the 
WKB approximation breaks down near the classical turning points. 
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factors inside the product (33) actually lead to a product of 6(= 2 + 1 + 1 + 2) 
factors, each of which is equal to a ratio of two T functions. Summary: we 
obtain a ratio of twelve T functions. 



The resulting condition conveniently written in terms of a shifted fre- 
quency / ~ 0(1) where —2if = —2iu — [—2ioj\ becomes 



' +M+2if+j ^-TY +-V/"+2v7'-y ^rv --^-/"+2v7-l \ 
±7 = ^ ^ ^ ^ 2 ' X f34^ 

-M+2if+j ^^^ -M+2if-j ^j^^ +M+2i/+l -J V-^-i; 



-.U+27/+J + 1 -^pj- -,'\.f+2//-,; + l -^pj- +,'\.f+2;/+2 ^ 
p^ +M+2i/+j+l ^p^ +M+2^/-i+l ^p^ -M+2^/+2 n 



(35) 



The factors n from (23) cancelled. Among these twelve F functions, six have 
argument with a huge negative real part and can be converted into T of 
positive numbers due to the well-known formula 

T(x) = — —-^^ (36) 

^ ^ sin(7rx)r(l -x). ^ ^ 

The factors of vr cancel, just like the Stirling approximations for the T func- 
tions with a huge positive argument: 

r(n + l) «i V27m(^)". (37) 

Only the sin(x) factors survive. The necessary condition for the regular 
asymptotic frequencies (i.e. frequencies where our analysis is valid; we will 
explain this point later) becomes 

^ sin7r(i/ + l) sin7r(i/ + |) sin7r(i/ - |) 

sin7r(i/+i) sin7r(i/+i±i) sin7r(i/ + i^) ^ ^ 

We chose M G 4Z so that we could erase M from the arguments of the 
trigonometric functions. We are also allowed to replace / by a; again because 
the functions in (38) are periodic with the right periodicity and we could have 
chosen the number [—2^0;] even. The terms (7r/2) in the denominator can 
be used to convert the sin functions into cos. The three arguments happen 
to be the same like in the numerator and we can also rewrite (38) as 

±i = tan(7riu;) tan(7riu; + 7rj/2) tan(7ria; — 7rj/2) (39) 

For generic complex values of j, (39) has formally six generic nontrivial 
solutions in every strip L < Im(aj) < L + 1 of the complex plane: three 
roots for each sign on the left hand side of (39). However two of them are 
infinite, l^j — > oo; the remaining four can be rewritten as exp(±47ra;) = 
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— 1 — 2cos(7rj).^ This form of our result can be easily derived if we multiply 
(38) by the denominator of the right hand side and expand the sin functions 
in terms of the exponentials. One can be more careful about the signs to see 
that the correct result must satisfy 

exp[e( Recu) ■ Attu] = -1 - 2 cos(7rj) (40) 

where e(y) is the sign function. Note that the modulus of the left hand side 
is never smaller than one, and therefore the asymptotic quasinormal modes 
do not exist for cos(7rj) < (—1/2). 

Let us now study (40) for various physical values of j. If j is a half- 
integer, the right-hand side of (40) equals (—1). The allowed frequencies are 
therefore^ 

jGZ + i: a; = ^^^^^^ + 0(n-V2), neZ (41) 
If j is an odd integer, the right-hand side of (40) equals (-1-1) and we obtain 
j G 2Z + 1 : ^ = y + 0(n-^/2), n G Z (42) 

Because the frequencies in (41) and (42) lead to an indeterminate form in 
(38) — in other words, co(n) and/or C2(n) vanishes for certain values of n — 
our analysis, strictly speaking, breaks down for these frequencies and more 
effort is needed to decide whether these solutions are physical. We will 
discuss these issues in the next subsection. At any rate, all these possible 
frequencies have a vanishing real part. 

On the other hand, the result for even j (scalar or gravitational per- 
turbations) is more interesting. There is no cancellation and (40) simplifies 
to 

JG2Z: a; = ^^^^^ ± + 0(n-n neZ (43) 

This is the solution that has been used to support loop quantum gravity. 

A partial summary: The only "regular" quasinormal frequencies that 
follow from our equations — without encountering indeterminate forms — are 
those with the real part equal to ± ln(3) (43) for j G 2Z. It is also conceivable 
that there are solutions with u = in/2 for odd j and u = i{n — l/2)/2 for 

*This form of the result was obtained in collaboration with Andrew Ncitzko. [19] 
®I am indebted to C. Herzog for pointing out an oversimplification in the first version 
of this paper. 
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half-integer j (with vanishing real part) where our analysis breaks down. 
In the next subsection we will show that these states indeed do exist. At 
any rate, =bln(3) is the only allowed nonzero asymptotic real part of the 
frequency. 

We cannot resist the temptation to say that the equation (40) simplifies 
not only for integer and half-integer values of j, but also for j E 2Z zt 2/3. 
In this (most likely unphysical) case we again obtain no solutions (even if 
one allows the sign function in (40) to be replaced by an arbitrary sign). 

3.2 Subtleties for purely imaginary frequencies 

The purpose of this subsection is to resolve the question marks regarding 
the special frequencies where our approximation, neglecting ci(ra), breaks 
down.^° 

First, we want to argue that the "regular" solutions are guaranteed to 
exist. Once we are able to relate the remainders R[-2iu)±M] through the 
continued fraction where we are allowed to neglect ci(n) as we said, we can 
also extrapolate them to (30). The boundary conditions Rq = oo, R^o = —1"*" 
imply that a specific sign of the square root in (30) must be separately chosen 
for n < [— 2iu;] and for n > [— 2ia;]. But the signs turn out to agree with the 
signs of ibi that lead to our solutions automatically. 

The necessary condition for u is also sufficient. The only solutions of 
this kind are the ln(3) solutions from (43). Their existence is guaranteed. 

In order to find all irregular solutions, we start with two useful observa- 
tions. 

• The continued fraction (28) depends on the coefficients co(n) and 
C2(n'), but only through their products co(n)c2(n -|- 1). 

• If zeroes appear in (33) exclusively in the numerator or the denomi- 
nator, the ratio R[-2iu>]-M / R[-2iu>]+M will be effectively either zero, or 
infinite. After we take ci(n) into the account, "zero" or "infinity" is 
replaced by a negative or positive power of respectively. 

The second point implies that the only possibility to obtain irregular 
solutions is to find oj such that (33) is an indeterminate form of type 0/0. 

'^"This subsection did not appear in the original version of the article. 
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Because co(n) can vanish at most for one value of n, there must be at least 
one value of n where C2{n) vanishes. Equation (23) then implies that 2iL0 + j 
or 2iu — j (or both) must be integer. 



For an integer or half-intcgcr spin, j, these two conditions are equivalent 
and therefore both numbers 2iuj it j must be integers for the quasinormal 
mode to have a chance to exist. In the half-integer case this condition pre- 
cisely agrees with the solutions proposed in (41): 

1 i(n- 1/2) 

iGZ + -: u; = ^^-^, neZ (44) 

Because the numbers 2iuj ± j differ by an odd integer for a half- integer j, 
one of the corresponding C2(n) factors appears in the numerator and the 
other appears in the denominator of (33) and allows for an indeterminate 
form. The coefficients co(n) arc then guaranteed to be nonzero for all n. 
The deviation of co from (44) affects the ratio in a generic fashion. For a 
suitable choice of this deviation, the indeterminate form is regularized to ±i 
once we include the effect of two relevant coefEcients ci (n) . We have checked 
the existence of solutions (44) with Mathematica. In the case of j = 1/2, for 
each allowed value of Im(u;) there are two solutions with nonzero =b Re(u;) 
that converge to zero. 



Now consider the case of integer j. As stated above (44), both num- 
bers 2iuj lb j must be integers. But if j is integer, then these two numbers 
differ by an even number. It follows that both vanishing factors of C2(n) 
appear in the numerator of (33) , or both appear in the denominator of this 
fraction. Without loss of generality, let us assume that they appear in the 
denominator. Then the indeterminate form is obtained only if the vanishing 
co(n) appears in the numerator. A closer look at (21), (23), and (33) implies 
that 2iuj lb j and 2iuj must be different modulo two, and therefore j must 
be odd. We can again find such a deviation from (42) that the effect of 
ci(n) gives us any desired result. For j G 2Z + 1 we therefore confirmed the 
existence of the states (42) while for j G 2Z the regular states from (43) are 
the only solutions. Our numerical calculations done using Mathematica fur- 
thermore indicate that for each allowed value of lm{oj) there are two states 
with nonzero values of =bRe(a;) that approach zero as Im(u;) — ^ oo. 



To amuse the reader, we would like to mention that for the special un- 
physical choice of the spin j G 2Z ib 2/3 discussed previously, at most one 
of the numbers 2iuj, 2iuj + j, 2iuj — j can be integer and therefore we ob- 
tain only one zero in (33) either in the numerator or in the denominator. 
This is not enough to end up with the indeterminate form of type 0/0. 
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Therefore for j G 2Z ib 2/3 there are no asymptotic quasinormal modes: 
neither regular nor irregular. This contrasts with the behavior for general 
j where we find four finite regular asymptotic quasinormal modes in every 
strip L < Im(a;) < L + 1. 

We showed that all irregular solutions that we have found are precisely 
the solutions where (39) led to a cancellation of zero against infinity. All 
solutions can be summarized by the equation (40) that we mentioned previ- 
ously. 

3.3 Generalizations and alternative computations 

It would be interesting to generalize our general procedures to related cases 
and apply the methods of answer analysis to other calculational frameworks. 
We have the following points in mind: 

• Reissner-Nordstr0m charged black holes. It is possible to write 
down the wave function as an expansion in {r — r^)/{r — r^) similar to 
(19) [12, 24]. The analysis of the critical behavior could be analogous 
to ours. One might be forced to add one or two terms to (20) because 
the Reissner-Nordstr0m warp factor contains an extra term {Q'^ /r"^). It 
has been recently shown by functional [19] and numerical [24] methods 
that the universal Hod-like result [8] 

Re{ui) ln(3)THawkmg as Im(a;) oo. (45) 
is replaced by a more complicated equation: 

/Hawking'^ + 2 + 3e''''^Hawking'^ = 0, k=— (46) 

r+ 

• Kerr and Kerr-Newman rotating (and charged) black holes. 

In the case of Kerr black holes (with Q = 0) , it is still possible to sep- 
arate the differential equations into radial and angular equations; see 
for example [12, 25]. Investigation of the quasinormal modes by the 
methods of this paper as well as the methods of [19] is in progress [26]. 
But some numerical calculations have already been done [24]. They 
indicate that the answer ln(3)rHawking is not universal; in fact, their 
calculation seem to imply that the quasinormal modes do not exist for 
a > acrit where the critical angular momentum is acrit ~ This 
result would mean that the connection between the area quantum and 
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quasinormal modes of the Schwarzschild black holes [6] was a coinci- 
dence that does not generalize to other black holes. Therefore it is 
desirable to confirm or reject the results of [24] by other means [26]. 

• Black holes in AdS spaces. A more detailed analysis of black hole 
quasinormal modes in Anti de Sitter space might be interesting espe- 
cially because the Maldacena's AdS/CFT correspondence might shed 
a new light on their origin. Papers [27, 28] show that the facts about 
the black hole quasinormal modes in the AdS space can be reinter- 
preted in terms of thermal physics of the dual Conformal Field Theory. 
In [29], the continued fractions play an important role. However, the 
asymptotic structure of the quasinormal modes of the AdS black holes 
differs from the flat space in many important respects. These differ- 
ences follow from a different asymptotic geometry of the spacetime. 
For example, even the asymptotic frequencies depend on I. Not only 
the imaginary part, but also the real part of the frequencies diverge. 
There is no ln(3) limit. The recurrence relation describing the quasi- 
normal modes of AdS black holes has also been studied [30, 31, 32]. 

• De Sitter space and Rindler space. The usual perturbations of 
de Sitter space form an integrable system that does not exhibit the 
ln(3) effect, at least in the simplest models that we checked (which 
however corresponded to j = 1 where our prediction for Re(a;) van- 
ishes anyway). The perturbations of Rindler space can be rewritten 
as a quantum mechanical model with the potential exp(x). This is 
the asymptotic behavior of (2) for x — > — oo and a general property 
associated with the horizons. The solutions are simply Bessel func- 
tions; see [33] for another application of the same Bessel equation. We 
do not know a clear argument determining what the boundary con- 
dition at X ^ +00 should be. Therefore we cannot derive a similar 
ln(3) effect in the Rindler space; it is even impossible to associate a 
unique temperature scale with the whole Rindler space and therefore 
it might be unreasonable to expect any definite answer. Another con- 
tradiction is that the result should be j-dependent as we have seen in 
the Schwarzschild case, but the Regge- Wheeler-like potential for the 
Rindler space seems to be ^-independent. 

• Higher-dimensional black holes. The calculation of the higher- 
dimensional Schwarzschild black hole's quasinormal modes could lead 
to much more important and interesting results. This is not just our 
belief; see the recent speculative attempt [34] to predict the higher- 
dimensional behavior. Even if the general structure of the modes re- 

^^I am grateful to V. Cardoso for reminding me of the Anti de Sitter space. 
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mains valid, the factor of ln(3) in (45) could be replaced by a different 
constant — for example ln((i — 1). The recent paper [35] showed that 
it was possible to calculate quasinormal modes in higher dimension; 
their asymptotic behavior was finally calculated in a very recent pa- 
per [19] that is based on functional methods. The quasinormal modes 
of scalar perturbations of Schwarzschild black holes seem to approach 
ln(3)rHawking in any dimension. This result seems puzzling from the 
viewpoint of loop quantum gravity [36]. It is desirable to reproduce 
or reject this result by numerical calculations and/or by the methods 
of this paper. However, the fivc-dimcnsional Schwarzschild black hole 
would probably lead to a recursion relation of the fifth order that is 
much more difficult to handle. 

Let us mention a very unlikely (but not quite impossible) speculative 

scenario. Imagine that for 11-dimensional black holes or black branes 
of some kind, ln(3) in (45) is replaced by ln(248) in the case of gravi- 
tational perturbations. Such a coincidence would probably encourage 
many to search for a spin network description of (the bosonic part of) 
the 11-dimensional supergravity based on the Es gauge theory. Such 
a description might be related to the usefulness of gauge theory in 
the bulk for the path integral quantization of the three-form potential, 
studied by Diaconescu, Moore, and Witten [37]. The Eg gauge field 
has a sufficient number of components to include the three-potential 
(expressed as the Chern-Simons three-form) as well as the metric that 
could be perhaps written in the (almost) usual LQG fashion: 

248 ^ r 

(det 510x10)5"' = -J2sA^s¥ «' ^ = 1' 2, . . . 10. (47) 

1=1 " b 

The field has many components, 10 X 248. While it might sound 
like a very redundant choice, we think that the meaningful proposals 
of LQG can be generalized to any spacetime dimension, as long as we 
allow the gauge theory configuration space to be bigger than the con- 
figuration space of pure gravity. For example, the quantization of two- 
dimensional areas in four spacetime dimension must generalize to the 
quantization of {d — 2)-dimensional areas in d spacetime dimensions — 
which is directly implied e.g. by (47). The reason is simply that the 
{d — 2)-dimensional areas determine the entropy. 

• The WKB approximation. The article [18] tried to solve the prob- 
lem of the highly damped quasinormal modes in the WKB approxima- 
tion. Although the authors could calculate the subleading corrections 
at very high orders, it was later realized that a subtlety invalidates this 
WKB approximation even at zeroth order because their asymptotic 
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Rc(a;)Guiim = for the highly damped modes is incorrect — certainly 
for general real values of j. Subsequently Andersson and Linnaeus 
[38, 14] improved the method of [18] and looked for the highly damped 
quasinormal modes again but the asymptotic behavior was not under- 
stood analytically. We should note that the WKB prescription of [18] 
seems to break down in a controllable way. Furthermore, the resulting 
frequencies should be a periodic function of j that appears, through 
the combination (l—j^), in the coefficient of the (1/r^) term of the po- 
tential in [18]. The correct result was finally reproduced by functional 
methods in a very recent paper [19]. 



4 Speculations on implications for gravity 

We start this section by rewriting the equation (8) for the quasinormal modes 
in usual units where 2M ^ 1 and Gn i^'^'- 

i(n - 1/2) ln(3) ^, _1/2^ 

Note that the ratio between the real part and the spacing of the imaginary 
part equals ln(3)/27r, regardless of the choice of units. It is useful to recall 
that the Hawking temperature of the Schwarzschild black hole of mass M 
equals Tnawking = l/(87rGiv-^)- Consequently, the equation (48) can be also 
written as 

— = 27ri(n-l/2) + ln3 + 0(ra-^/2-) ,^49^ 

Hawking 

We can remove the term proportional to n simply by exponentiating (ttz in 
(49) gives the minus sign below): 

expf-^^) =-3 + 0(n-V2) (50) 

\ Hawking / 



Because the quasinormal modes show the position of poles of the transmis- 
sion amplitude in the corresponding quantTim mechanical model — and these 
poles will probably also have some interpretation in the spacetime terms — 
the asymptotic structure of this transmission amplitude (for large Im(a;), 
and perhaps for large |a;| in general) will be proportional to the factor 

^('^) ~ Att ^ uj ^' /^Hawking = Tf^-^ • (51) 

gPHawkmg'^ ^ 3 i Hawking 



The roots obtained by the reflection symmetry have a negative imaginary 
part — ln(3) and lead to a similar factor l/(cxp(/^Hawkingt^) + 1/3)- These 
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denominators are reminiscent of thermal physics. Actually, had we started 
with our results for the odd spin (42) and the half-integer spin (44) , the same 
procedure would have given us a more familiar factor in the transmission 
amplitude proportional to 

'^^^^ ^ e%awking'^ :p 1 ■ ^^^^ 

This is exactly the average occupation number in Bose-Einstein statistics (for 

j odd the denominator has a minus sign) or Fermi-Dirac statistics (for j half- 
integer the denominator has a plus sign). In both cases the statistics agrees 
with the spin j of the perturbation. We can imagine that the transmission 
amplitude results from an interaction of the given particle with a thermal 
bath of temperature Tkawking) containing particles of the same statistics. In 
other words, the amplitude is related to the thermal Green's function in the 
Schwarzschild background; this explains the general form of (51) and (52). 

This agreement makes the result (51) for the even values of j even more 
puzzling. Why do we fail to obtain the same Bose-Einstein factor as we did 
for odd j7 Instead, we calculated a result more similar to the half-integer 
case, i.e. Fermi-Dirac statistics with the number 3 replacing the usual number 
1; let us call it Tripled PauU statistics. Such an occupation number (51) can 
be derived for objects that satisfy the Pauli's principle, but if such an object 
does appear (only one of them can be present in a given state), it can appear 
in three different forms. Does it mean that scalar quanta and gravitons near 
the black hole become (or interact with) J = 1 links (triplets) in a spin 
network that happen to follow the Pauli's principle? Our puzzling results 
are summarized in the table below. 



Spin 


Asymptotic frequencies 


Corrcspontling poles 


Naively implied statistics 


j G Z + 1/2 


i(n-l/2) 
2 


1 

*llawldii2;-') + l 


Fcrmi-Dirac 


i G 2Z + 1 


m 
2 


i 

exP(/3HawkinK'^)-l 


Bose-Einstein 


jG2Z 


± ln(3) 


1 

exp(/3HawkinK'^)+3 


Tripled Pauli? 



(53) 



4.1 Path integrals and speculations on the black hole chem- 
ical potential and spin network resonances 

It is a well-known fact that the asymptotic periodicity of the time coordinate 
in the Euclidean black hole solutions equals the inverse Hawking temperature 
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/^Hawking exactly if there is no deficit (or excess) angle at the horizon (which 
is a "tip" of the solution). In this sense, the Euclidean solution "knows" 
about the Hawking temperature. Black holes can only be in equilibrium 
with a heat bath of the Hawking temperature because if we want to consider 
black hole configurations contributing to the path integral, we are forced to 
take the right periodicity of the time direction. 

If we glance at (52), we see that the quasinormal modes also "know" 
about the right temperature. The asymptotic quasinormal frequencies in 
the case of spin j ^ 2Z satisfy exp(/3Hawking'^) = ±1- But exp(/3Hawking'^) 
is the evolution operator by /^Hawking in the Euclidean time. So the al- 
lowed asymptotic frequencies are exactly the frequencies that respect the 
time periodicity. Alternatively, the transmission amplitude (52) contains 
the information about the thermal ensemble. Both Hawking radiation as 
well as quasinormal modes deal with general relativity linearized around the 
black hole solution, so it is perhaps not too surprising that they contain a 
similar piece of information. 

We understand that (52) is a natural expression related to the thermal 
properties of the black hole. But what about (51), which has +3 instead of 
— 1 in the denominator? These quasinormal frequencies appear for j G 2Z. 
They satisfy exp(/?Hawking'^) = —3 (or —1/3, if we consider the solutions 
with the negative real part). Just like the solutions (42) and (44) lead 
to (52) which correspond to the path integral with periodic or antiperiodic 
boundary conditions on the time coordinate of periodicity /^Hawking; the extra 
frequencies (43) seem to lead to the path integral where the fluctuations 
get multiplied by (—3) (or (—1/3), which is related to (—3) by the time 
reversal symmetry) if we perform a rotation around the Euclidean time. Such 
boundary conditions correspond to the evaluation of the thermal expectation 
value of the operator 

((-3)^) = IV [exp(-/3Hawking^)(-3)^] (54) 

where N is the number of gravitons (or scalar quanta) in the state. Just 
as the frequencies (41) and (42) inform us that the thermal expectation 

values of various ordinary operators contain the denominator (52), the extra 

quasinormal modes (43) might perhaps hide a similar insight. Note that the 

insertion of (—3)^ in (54) might be interpreted as a chemical potential^^ for 

gravitons (or scalar quanta) in a grand canonical ensemble e^p{—(3H + iJ,N). 

Something special should happen when the chemical potential /x formally 

approaches the (complex) value 27ri(n + 1/2) + ln(3). 

^^We are grateful to Alex Maloney for suggesting the possible relevance of the term 
"chemical potential" in this context. 
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In loop quantum gravity, In (3) arises as the entropy of a single link 
with the minimal possible spin 

j(min) ^ I 

assume the S0{3) version of 
LQG). An advocate of loop quantum gravity could argue in the following 
way: (—3)^ is an operator that receives some "resonant" contributions from 
the nontrivial quasinormal modes. The reason could be that (—3)^ might 
actually be the operator that creates (or destroys) a single link of the spin 
network. The number 3 (or (—3), we have no understanding of the minus sign 
at this point; it might arise from Fermi-like statistics of some objects such as 
the J = 1 links) would be related to the number of possible Jj polarizations 
of the new link. The number of gravitons N would be correlated with the 
number of the links in some unknown way. Finally, there could also be 
an explanation why the creation operator of a link has a large resonant 
expectation value in the canonical ensemble describing a black hole. Only 
time can show whether these speculations can be supported by a consistent 
formalism. 



5 Conclusions and outlook 



In this paper we calculated the asymptotic form of the frequencies of the 
quasinormal modes on the four-dinicnsional Schwarzschild background. We 
used the method of Taylor expansion and continued fractions that can be 
approximated very well for large frequencies by imposing the constraint that 
the quantities Rn, describing the continued fractions, change slowly. For the 
values of index n near \n + 2iuj\ ~ 1, this estimate breaks down for general 
values of lj. But for large |u;| the recursion relations near \n + 2iu;\ « 1 can 
be approximated by an exactly solvable system. The continued fractions 
become ordinary fractions. The expressions can be rewritten in terms of the 
r function and then even in terms of trigonometric functions. The result for 
odd and half-integer spins j of the perturbation are the frequencies 

= ^Hawking • 2m{k + j), keZ (55) 

but for even integer spins j we proved the existence of modes with 

= rHawking[27ri(fc + 1/2) ± In 3], k e Z. (56) 

We conjectured that these asymptotic frequencies of quasinormal modes can 
be found for other black holes, too, as long as the frequency is written as a 
multiple of the Hawking temperature; the result (56) was recently confirmed 
in the case of higher-dimensional Schwarzschild black holes [19], but ruled 
out in the Reissner-Nordstr0m case [19] as well as the case of Kerr black holes 
[24]. While (55) carries some information on the thermal structure of the 
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black hole, the significance of the nontrivial frequencies (56) remains clouded 
in mystery. We presented speculations how (56) might be understood in the 
language related to loop quantum gravity. 

It would be interesting to 

• check more carefully whether all a priori allowed solutions in (55) and 

(56) exist; look for the "irregular" solutions (55) numerically; it is 
however true that some very difficult numerical calculations become 
less interesting because their computational goals can now be achieved 
analytically 

• generalize the calculation to higher-dimensional black holes and derive 
what ln(3)THawking is replaced by; a very recent paper [19] confirmed 
that ln(3)rHawking appears for scalar perturbations of Schwarzschild 
black holes in an arbitrary dimension 

• try to apply similar methods to the case of charged black holes and 
black holes in the Anti de Sitter space; functional methods have been 
used to derive the aperiodic structure of the Reissner-Nordstr0m black 
holes in four dimensions [19] and this structure has been confirmed in 
a very recent paper [24] 

• investigate the asymptotic structure of the quasinormal modes of ro- 
tating Kerr black holes in four dimensions; preliminary results suggest 
that ln(3)rHawking is not the universal answer in this case [24] 

• calculate the sublcading contributions in l/( Im(aj))^/^ to Rc(lj) using 
our methods or the methods of [19] as the zeroth order approximation 

• find an interpretation of the real part that does not involve loop quan- 
tum gravity 

• or even more ambitiously, find a solid interpretation in terms of loop 
quantum gravity; our analysis never led to the number ln(2); does it 
show that loop quantum gravity is inconsistent with the existence of 
fermions? 

Obviously, there are still many tasks to be solved. 
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